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Abstract 

This  report  presents  a  proof  of  the  uniqueness  of  a 
parallel  three-dimensional  shear  flow  in  a  channel  with 
arbitrary  cross  section  where  the  speed  of  the  flow  Is  not 
less  than  the  highest  critical  speed.   The  Investigation 
also  Includes  a  two-dimensional  analysis  In  which  It  Is 
assumed  that  while  the  flow  velocity  varies  with  the  depth, 
the  density  also  depends  on  the  depth;  and  for  this  case 
the  development  leads  to  a  formula  which  gives  a  good 
approximation  to  the  highest  critical  speed. 
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1.   Introduction 

Consider  an  Incompressible,  Invlscid  liquid  contained 
In  a  horizontal.  Infinitely  long  straight  channel  whose  cross 
section  Is  arbitrary.   Suppose  that  a  gravitational  force  Is 
the  only  body  force  which  acts  on  the  liquid.   The  nonlinear 
hydrodynamlcal  equations  given  below  In  Section  2  show  that  a 
uniform  parallel  flow  Is  a  possible  steady  motion.   This  kind 
of  flow  Is  defined  to  be  such  that  the  only  non-zero  velocity 
component  Is  the  component  v.  In  the  axial  direction  of  the 
channel;  and  v-.  ,  although  assumed  to  be  Independent  of  the 
longitudinal  coordinate,  may  be  a  function  of  the  lateral  co- 
ordinates of  the  channel.   If  the  equations  are  linearized  with 
respect  to  a  certain  parallel  flow  the  resulting  linear  equations 
also  admit  a  similar  flow  and  in  particular  the  uniform  parallel 
flow  in  which  the  axial  velocity  is  constant.   However,  according 
to  the  linear  theory,  this  is  not  the  only  possible  motion  if 
the  speed  of  the  liquid  at  infinity  is  less  than  one  of  a  possible 
set  of  critical  values.   For  example,  if  the  cross  section  of 
the  channel  is  a  rectangle  with  depth  h  and  if  the  speed  of  the 
liquid  at  infinity  is  less  than  the  critical  speed  /gh,  where  g 
is  the  acceleration  due  to  gravity,  then  the  linear  equations 
predict  that  a  progressing  wave  motion  is  possible. 

A  discussion  of  critical  speeds  is  necessary  for  the 
analysis  of  several  hydrodynamlcal  problems  concerned  with 
channel  flow.   They  arise  in  the  study  of  the  motion  due  to  a 


surface  pressure  distrubance  which  moves  In  the  direction  of 
the  channel  with  fixed  speed  either  when  this  problem  is 
regarded  as  a  steady  state  problem  or  when  it  is  regarded  as  a 
Newtonian  initial  value  problem.   In  the  steady  state  analysis 
of  the  problem  critical  speeds  arise  not  only  with  respect  to 
the  uniqueness  of  the  solution  but  also  with  respect  to  the 
admissibility  of  the  linearization.   In  the  Newtonian  approach 
based  on  an  Initial  value  problem  for  the  linear  theory  it 
turns  out,  as  Stoker  [1]  showed,  that  at  a  critical  spped  the 
velocity  components  of  the  flow  become  unbounded  as  time  elapses. 
The  nonlinear  theory  of  a  gravitating  fluid  in  a  channel  leads 
to  the  interpretation  of  critical  speeds  as  bifurcation  values 
at  which  cnoidal  and  solitary  waves  may  appear  as  well  as 
parallel  flows.   These  examples  point  to  the  fact  that  critical 
speeds  can  be  defined  in  different  ways.   A  discussion  of  the 
various  definitions  can  be  found  in  a  paper  by  Peters  and 
Stoker  [2]. 

During  conversations  with  the  author  and  about  problems 
similar  to  those  mentioned  above  J.  J.  Stoker  raised  the 
following  uniqueness  question.   If  a  gravitating  liquid  confined 
to  a  rectangular  channel  is  in  a  state  of  parallel  flow  with  a 
finite  speed  not  less  than  the  highest  critical  speed  does  the 
linear  theory  show  that  this  flow  is  the  only  possible  steady 
motion  which  is  bounded?   In  the  sections  which  follow  we  show 
that  the  answer  to  this  question  is  in  the  affirmative.   We  show 
this  under  the  assumption  that  the  density  of  the  liquid  varies 


with  depth  and  that  the  liquid  is  subject  to  a  shear  In  velocity. 
Our  method  is  based  of  course  on  an  eigenvalue  problem  which 
possesses  only  the  trivial  solution  provided  that  a  parameter  of 
the  problem  is  not  less  than  a  certain  value. 

Weinstein  [5]  showed  that  If  (j)(x,y)  is  a  potential  function 
which  is  required  to  satisfy 


!•  ^xx^^'^^  "^  ^yy   =  0  ,         -oo  <  x  <  oo 


0  <  y  <  1 


2.  (tytx^O)  =  0  , 


5.  <t>y(x,l)  =  P<l)(x,l)  ,  p  >  0 


and  if  A   is  the  unique  positive  root  of 


A   tanh  A  =  p  , 
o       o   '^  ' 


where  p  is  a  constant,  then 


Mx.y)  =  [a   cos  A  +b   sin  A  xlcosh  A  y 

T  \      'J  '  Lq         oo        O  O 


is  the  only  bounded  function  which  satisfies  the  above  conditions, 
Weinstein 's  proof  of  this  is  based  on  a  completeness  theorem. 
In  Section  3  of  this  paper  we  use  Weinstein 's  method  to  analyze 
the  eigenvalue  problem  which  we  derive  in  order  to  discuss  the 
two-dimensional  flow  of  a  gravitating  liquid  with  non-constant 


density  and  velocity  each  varying  with  depth.   In  the  course  of 
the  analysis  we  find  a  formula  which  gives  an  approximation  to 
the  highest  critical  speed. 

Section  h   Is  devoted  to  an  analysis  of  a  three-dimensional 
motion  of  a  gravitating  liquid  of  constant  density  In  which  the 
velocity  depends  on  the  coordinates  orthogonal  to  the  direction 
of  the  containing  channel  which  Is  assumed  to  have  an  arbitrary 
cross  section.   The  character  of  the  eigenvalue  problem  which  we 
formulate  for  this  case  Is  different  from  that  presented  for  the 
two-dimensional  case.   As  a  consequence.  Instead  of  seeking  a 
method  based  on  a  completeness  theorem,  we  base  the  analysis  on 
the  generalized  Fourier  transform  theorem  which  incidentally  can 
also  be  used  for  the  case  of  Section  5  in  lieu  of  Welnsteln's 
procedure. 


2.   Formulation 

Let  a  gravitating,  incompressible,  and  inviscid  liquid  with 

density  p  be  confined  to  an  infinitely  long  horizontal  channel 

whose  cross  section  is  constant.   Suppose  that  the  equilibrium 

free  surface  of  the  liquid  is  planar  and  that  it  coincides  with 

the  horizontal  x-,  ,x-,-plane  of  a  cartesian  reference  frame  whose 

1  y 

x-,-axls  is  taken  parallel  to  the  rigid  cylinder  which  forms  the 
channel.   With  the  positive  direction  of  the  Xp-axls  taken  to  be 
upward,  let  the  channel  wall  be  defined  by 


x^  =  Q(x^) 


and  let  the  free  surface  be  given  by 


Xg  =  F(x^,x^,t)  . 


Let  g  denote  the  gravitational  acceleration,  let  ir  denote  the 
pressure,  and  let  us  use  v, ,  v„,  v   to  denote  the  velocity 
components  of  a  liquid  particle,  while  t  stands  for  time.   In 
terms  of  these  quantities  the  elementary  theory  of  hydrodynamics 
predicts  that  if  the  gravitational  force  is  the  only  force 
acting  then  the  motion  of  the  liquid  is  defined  by  the  continuity 
equation 

Sv    hv  bv 

(2.1)  ^  +  ^  +  ^  =  0  , 

dx-L    0X2    ox^ 

the  incompressibility  condition 

(2.2)  V,  |£-+  v„  |£-+  v^|£-=  0  , 

1  ^x^     2  ^xg     3  dx^ 

the  momentum  equations 

^^^t  1   ax^  2   dx^  3   dxT  hx-^    ' 

Sv  Sv  Sv  av  > 

(2.3)  p(^  +   V,         2   ^  2   ^  2)    _   pg   _   dTT 


1 

hx-^ 

"2 

hx2 

1 

dx^ 

+ 

^2 

3vj 

Sxg 

^^  -   '1  ^^  -  '2  ^^  -  '3  -^^   -  -  ""^  -  ^ 


Sv-,  Sv-.  hv-,  hv-,  -^^ 

^'at         1  dxZ       2  axo       3  hxJ        ^xZ 


the  kinematic  boundary  conditions 


i       :5 
the  dynamic  free  surface  condition 

(2.6)  7r(x^,F,x^,t)  =  0  , 

plus  initial  conditions  at  t  =  0,  and  conditions  which  specify 
the  behavior  of  the  liquid  at  distances  arbitrarily  far  from 
the  origin. 

The  above  equations  can  be  written  in  dimensionless  form 
if  we  introduce  a  typical  length  in  the  vertical  direction,  say  h, 
and  the  dimensionless  quantities 

v,-l  u-1  K-1 

X  =  x-j^h        ,    y  =  x^h        ,     z  =  x  h        , 

^1  "  ^i^gh)~     '    ^2  "  ^2^S^^~     '        ^3  "  v^(gh)"  /  , 

TT^  =  7r(pgh)"^    ,    f^  =  Fh~^       ,     q  =  Qh"^ 

T  =  t(g/h)l/2   ^ 

where  p  is  some  fixed  quantity  with  the  dimensions  of  density. 
In  terms  of  these  quantities  the  equation  of  the  channel 
wall  is 


y  =  q(z)  , 

the  equation  of  the  free  surface  Is 

y  =  f^(x,z,T)  , 

and  the  basic  hydrodynamlcal  equations  are 


Su-,    ^Up    ^u^ 


Su,      Su-.      ^u,      Su^        Stt, 

SUp        SUp        SUp        5Up  ^TT, 

(2.9)     p(^  +  u^^  +  U2^+u^^)  =  -  p-^^  , 

Su^      Su--      Su^      Su         Stt^ 
•^Mt     1  dx     2  dy     3  dz       ^  dz   ' 


with  the  boundary  conditions 


(2.10)  u^  =  u 


Sq 
2  =  "5  ^ 

"br  Sf   Sf 

(2.11)  U„  =  U^  -r-=-   +  U-,  -^-i  +  3-^ 

^     '  2    1  dx     3  dz     dT 


and 


(2.12)  •n-^(x,f^,z,T)  =  0  . 


} 


This  system  is  satisfied  by  the  quantities 


u^  =   u^(y,z)  = 


7  +VQ(y,z) 


(2.13)     /  u^  =  0  ,    u^  =  0  ,    f-^  =  0  , 


p  =  ppQ^y)  '     ^1  "j    pQ^^^'^^  ' 


where  /  is  constant  and  v   is  a  continuous  non-negative  function. 
They  define  a  steady  parallel  motion  in  the  channel  and  we  will 
refer  to  this  flow  as  the  equilibrium  flow.   The  velocity  v  {y,z) 
gives  the  transverse  shear  in  the  axial  velocity  component  and  it 
is  also  a  measure  of  the  departure  of  the  flow  from  a  uniform 
state  defined  by  the  velocity  7.   The  function  pAv)   measures  the 
variation  in  density  with  the  depth  and  we  suppose  that  it  does 
not  decrease  as  the  depth  Increases  so  that  with  respect  to  our 
coordinate  system  the  derivative  p'{y)  if  it  exists  satisfies 


dPo(y) 


dy   - 


<  0 


Let  us  proceed  to  linearize  the  equations  (2.7)-(2.12)  with 
respect  to  the  flow  given  by  (2.13).   That  is,  let  us  write 
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u^  =  u^(y,z) +  u  ,  ^  =  pQ^y)  + '^   > 


(2.14) 


U2  =  V  , 


0 


u^  =  w  , 


f^  =  f(x,z)  , 


and  assume  steady  motion.   Let  us  substitute  these  quantities  In 
equations  (2.7)-{2.12)  and  neglect  terms  which  Involve  products 
of  two  or  more  factors  from  the  set  u,  v,  w,  0",  p  and  f.   The 
result  of  the  linearization  of  the  equations  (2. 7) -(2. 9)  Is 


(2.15) 
(2.16) 


u  +  v  +  w  =0 
X   y   z 


u    +  vp   =0 
o  X   "^oy 


(2.17) 


'^  p  fu  u  +U   V  +u  w)  =  -  p 
^o*  o  X   oy    oz  '     ^x 


<  p  u  V 
\  "^o  o  X 


p  u  w 

^O  O  X 


=  -  5--  p. 


-  p. 


The  condition  at  the  channel  wall  Is 


(2.18) 


V  =  w 


Sq 
^ 


With  respect  to  the  free  surface  conditions  (2.11),  (2.12)  they 
become  conditions  to  be  satisfied  at  y  =  0.  In  place  of  (2,11) 
we  have 


(2.19)  v(x,0,z)  =  u^(0,z)f^(x,z)  , 


and  from  (2.12)  we  have 

0 


J      p^{r])a^   +   p(x,f,z)  =  0 


which  after  differentiation  and  removal  of  second  order  terms 
becomes 


(2.20)  -p^(0)f^(x,z)  +  p^(x,0,z)  =   0 


Notice  that  if  v  =  0,  w  =  0,  the  linearized  equations  are  again 
satisfied  by  a  flow  of  the  type  (2.13). 

Our  object  now  is  to  show  that  if  the  speed  7  is  not  less 
than  a  certain  highest  critical  value  then  the  only  possible 
bounded  solution  of  the  problem  formulated  by  the  equations  (2.15)- 
(2.20)  is  the  one  which  defines  an  equilibrium  flow,  (2.13). 


3.   Two-dimensional  Motion.   Rectangular  Channel 

If 

7  +  ^^(y) 

^o  =  ^o^y^  =  - 


gh 


where  v  is  continuous;  if 
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Po  =  Po^y^  ' 


w  =  0 


and  If  the  remaining  quantities  in  the  equations  (2. 15)-(2. 20) 
are  independent  of  z,  then  these  equations  define  a  two- 
dimensional  motion  which  may  be  interpreted  as  a  two-dimensional 
flow  in  a  rectangular  channel.   For  this  case  the  basic  linearized 
equations  are 

(3.1)  u^+Vy  =  0 

(5.2)  u^cr^+vp^y^O 

(3«3)  P^C^^u  +u  v)  =  -  p 

*  f^o*  o  X   oy  '     '^x 

^o  o  X  y 

If  the  depth  of  the  channel  is  h  the  equation  of  the  bottom  in 
our  dlmensionless  variables  is  y  =  -1  and  since  the  vertical 
velocity  component  must  vanish  there  we  must  have 

(3.5)  v{x,-l)  =  0  . 

Corresponding  to  the  free  surface 

(3.6)  y  =  f(x) 
the  linearized  free  surface  conditions 
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(5.7)  v(x,0)  =  u^(0)f^{x) 


(3.8)  -Po^^^^x  +  Px^^'°^  "  0 


must  be  satisfied. 

For  the  analysis  of  the  two-dimensional  equations  we  will 
work  with  the  dependent  variable  v/u  which  we  denote  by  X(x,y) 
The  elimination  of  u,  p  and  CJ  from  these  equations  shows  that 


must  satisfy 


^^•9)      1^  (Po^O^X^  +  ^  (Po^oXy)  -  Poyt  =  0  ,      _1  <  y  _  0  , 


-OO  <  X  "^  00 


with  the  boundary  conditions 

(3.10)  7((x,-l)  =  0 


(3.11)  u^(0)  ty{^,0)  -  Xix,0)    =0 


We  turn  now  to  the  method  of  Weinstein  [3]  and  replace 


^  ^Po^o^x^ 


2    2 

in  equation  (3-9)  with  -A  p  u  fiy);    while  we  replace  '^   and  X 

respectively  with  f     and  f.      This  formulates  and  introduces  the 
following  standard  eigenvalue  problem: 


12 


(3.12)    ^  (Po^o^)  -Poy^fy^  =  ^^Po^o^  '  -1  <  y  <  0 


(3.15)  ^(-1)  =  0 


(3.1^)  u^(0)^y{0)  -^(0)  ^   0  . 


Here,  the  eigenvalues  depend  on  the  magnitude  of  y   and  the  elgen- 
functlons  must  satisfy  the  orthogonality  relation 


0     m  ^  n 

pQ(y)uQ(y)^ni(y)^j^(y)dy  = 

-1  1  1     m  =  n 


/ 


where  f     and  ^  are  the  elgenfunctlons  which  correspond 

respectively  to  the  eigenvalues  A  and  A  .   It  Is  easy  to  see 

n      m 

that  If  A   Is  an  eigenvalue  then  so  Is  -A  .   Besides  this.  It  Is 
n         °  n 

not  difficult  to  verify  that  the  eigenvalues  A  are  either  real 
or  pure  Imaginary  numbers.   It  Is  known  that  there  can  be  only  a 
finite  number  of  eigenvalues  In  any  bounded  domain  of  the  complex 
A-plane.   It  Is  also  well  known  that  the  above  second  order 
system  defines  a  complete  set  of  elgenfunctlons  {^„(y)}  such  that 
a  twice  dlfferentlable  function  G(y)  can  be  expanded  In  the 
Fourier  series 

00 

G(y)  =  XZ  "n^n^y)  • 
n=0 
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It  follows  from  the  above  remarks  that  for  any  fixed  value 
of  X,  ')(l(x,y)  has  the  unique  expansion 

X   (x,y)  :=  ^^   =  ZZ  ^n^^)^n^y)    >  "1  1  Y  1  0  • 

o^*^  "^    n=0 


Here,  the  coefficient  ct^fx)  is 

0 

a^{x)  =J    PQ(y)u^{y)X(x,y)^n^y^^y 
-1 

and  as  a  function  of  x  it  must  satisfy 


dx      _-, 


V-  =  j  po^o^xx^y 

0 

=  -/  [|rPo-oXy-Poy^l^n^y 
-1 

0 

=  -  ^n  /  Po^o-X  V^  • 


-1 


Hence  the  coefficient  a^(x)  must  satisfy 


^ —  +  -K'^a    (x  =  0 

,  2  n  n^  ^ 

dx 

from  which 

a^(x)  =  a^  cos  A^x  +  b^  sin  A^(x)  ,        \  ^  ^ 
(3.15) 


Ik 


This  shows  that  if  'X  (x,y)  Is  to  be  bounded  everywhere  then  If 

a.  does  not  correspond  to  a  real  eigenvalue  we  must  take  a  =  0. 
J  J 

If  A   =  0  Is  an  eigenvalue  the  corresponding  term  in  the  develop- 
ment oT    %    =   v/u  Is 


Since 


(^o+V^^o^y^  • 


U   =  -V 

X     y 


the  corresponding  term  In  the  development  of  u  Is 

2 

-(a  X  +b  4^-)^   +f^(y) 
^  o    o  ^  '  oy   o'*^ ' 


but  this  Is  unbounded  as  x  — >  od  and  therefore  for  a  bounded 
motion  we  must  take  a  =  0.   It  Is  now  apparent  that  our  problem 
Is  reduced  to  a  study  of  how  7  affects  the  disposition  of  the 
eigenvalues  of  the  system  of  equations  (3. 12)-(3. l4) . 

In  order  to  study  these  eigenvalues  let  us  convert  the 
equation  (^•12)  namely 


-^  P  u%   (y)  -  p  Tp  =   A^p  u^Tp   ,  -1  <  y  <  0 

dy  "^o  o  y      oy      ^00 


into  an  integral  equation.   If  we  integrate  (3. 12)  from  y  to  zero 
we  find 


15 


p^(o)u^(o)^y(o)  -PQ(y)u^(y)V/y(y) 


0 


0 


y 


+  A^  /    p^(Tl)u^(Tl)^(Tl)dTl 


from  which,  by  using  u  (O)'^  (O)  -f[0)    =  0,    we  obtain 

0 

PQ(y)u^(y)^y(y)  =  pQ{y)^(y)  +J     Po(ii)V'^(T])dTi 

y 

0 


-  A 


J       pQ^^^^o^^^^^^^^^  • 


y 


Since  ^(-1)  =  0,  the  last  equation  can  be  written 


y         0 
Po{y)uQ(y)^y(y)  =  p^Cy)  j   ^^{Ti)dTi+J  p^{r\)f^{r])dr] 


-1 


0  11 

-   A^y     Pq{ti)u^(ti)  jT    t^{|)dCdTi 


y 


-1 


which  after  an  integration  by  parts  yields 


y 


0 


Po(y)u^(y)^y(y)  =  Po(y)  /  ^  (Ti)dTi  +  /  p^{ti)^  (Ti)dT] 


-1 


y 


r 


0 


0 


-  A 


■^ 


p^(e)uj{^)d^  •  /    f^{i)di 


y 


y 


f  ir])  f   p^i^)ulii)dia^ 


0 
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or,    after  a  rearrangement. 


-1 

0  y 


2    I  y  -1 


-    A 

0  0 


+  /    \^^)'f     Poi^)^o^iW<ir\ 


n 
y  Ti 


If  we   set 

<t>  =  /po(y)  %(y)^y(y) 

the  last  equation  becomes 

/p^TyJ^iy)         {i  %{Ti)/p7rnT 

0  0 

^ f      ^^'^\  f      p      r^)u2(^)d^dTl 


+ 


Now  if  we  introduce  the  symmetric  kernels 


17 


k-,^(y,Ti;7)    = 


-1    <   T]    <   y 


y    <   r\    <    0 


0 


y 


UQ(y)uQ(Ti)/p^TyTpQ(Ti) 


ko(y,Ti;7)    = 


0 


n 

Uq  ( y )  u^  ( Ti ) /p^TyTprj^ 


-1   <   Ti    <   y 


y    <   T)    <    0 


we   find  the   Integral   equation 


0 


(5.16)        My)    =J      k^(y,Ti;7)<t)(Ti)dTi    -   ^    J      ^2^y,r\;7)H'^) 


-1 


0 
-1 


IdTi 


as  one  which  Is  equivalent  to  the  differential  system  (3.12) 

There  is  no  loss  of  generality  In  assuming  that 

0 


/  <t'^(y)dy  =  1  . 


-1 

Then  If  we  multiply  (3.16)  by  (f)(y)  and  Integrate  from  -1  to  0 
we  have 
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0       0 

r 

-1       -1 


1  +  A^J   My)/  k2{y,Ti)(t){Tl)dTidy  =  j  ^{y)  J    k^(y,Ti)(i)(Ti)dTidy  . 


P  0 

<t>(y)/ 
-1    -1 


The  Integral  on  the  left  hand  side  of  the  last  equation  is  equal 
to 


r   y 


-1 


,2 


f        Mn)dn 


_i  UQ(n)/p^TnT 


dy 


so   that 


0 

r 
-1 


r   y 


t2 


-1  %(^)/p^^TnT 

0  0 


dy 


=  /    <t>(y)  7    k^(y,'n)<l>(n)dTidy     . 


-1  -1 


An  application  of  Schwartz's  inequality  gives 


0  r     3 

l+A^jT    p^(y)u2(y)     y^ 


-1 


MTl)dTl 


.1   %ir\)/pj^ 


dy 


0        0 


-  /  /   /  ^i^y'^^^^^y 


-1   -1 


The  integral  of  k,  (y,Ti)  is 


0 


Po^y) 


-^IL dy 


-1  ^o^y^   -1  Po^^J^o^^^ 


Hence  we  see  that 
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(5.17) 


0 


r  y 


-1 


M'n)dTi 


-i2 


_i  Uq(ti)/p^TtT 


dy 


0   /  \   y 
Po(y)     r 


dT] 


-1  %^y^   -1  Po^^^^o^^) 


dy  -  1 


The  Inequality  (3- XT)  shows  that  there  can  be  only  a  finite 
number  of  real  eigenvalues  because  It  requires  these  to  lie  In  a 
bounded  line  segment  which,  as  we  noted  above,  cannot  contain  an 
infinite  number  of  eigenvalues.   The  inequality  (3-17)  also  shows 
that  if  7  in 

7+  v^(y) 


%M  ' 


O 


^ 


is  taken  sufficiently  large  then  the  absolute  magnitude  can  be 
made  as  small  as  we  please.  In  addition,  (3.17)  shows  that  if 
7  is  so  large  that 


0 


(3.18) 


Po(y) 


y 


dr) 


_i  %(y)    _i  Po(^)%(^) 


dy 


then  the  eigenvalues  cannot  be  real  and  non-zero.   From  what  has 
been  noted  above,  this  means  that  when  {3.18)  holds  the  only 
bounded  solution  of  {3.9)-(3.11)  is  the  trivial  one 
X^^^y)  ~   v(x,y)/u  (y)  =  0.   With  the  vertical  velocity  component 
V  equal  to  zero  everywhere  the  only  solution  of  the  original 
equations  {3.l)-(3.8)  for  the  two-dimensional  flow  is  one  of  the 
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type  (2.13).   In  other  words,  the  linear  theory  Implies  that  the 
flow 

7  +VQ(y) 

U-,  =  ,     Up  =  0  , 

ygh 

f^  =  0  , 

0 

p  =  ppQ^y)       >      "^1  "J    p(Ti)dTi 

y 

is  a  unique  steady  two-dimensional  flow  If  (5.18)  holds. 

In  the  foregoing  we  have  made  the  tacit  assumption,  which 
we  retain,  that  If  m  and  M  are  respectively  the  minimum  and 
maximum  values  of  the  non-negative  continuous  function  v  (y)  then 
y   does  not  have  a  value  between  -M  and  -m.   Without  such  an 
assumption  the  Integral 

y 

r  — *i^ 


-1  Po'l'"o'l' 


would  fall  to  exist. 

We  are  now  In  a  position  to  define  a  critical  speed  as  a 
speed  7  =  c   which  corresponds  to  a  transition  from  real  values 
of  A  to  pure  Imaginary  values.   This  occurs  when  A  passes  to  the 
zero  value  and  as  we  can  see  from  the  Integral  equation  {5.16) 
this  takes  place  when  7  =  c   is  an  eigenvalue  of  the  equation 

0 

(3.19)  My)  =  /   k^(y,Ti;7)M^)dTi  . 

"-1 
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It  can  be  seen  from  (3- 15)  that  c   Is  the  limit  speed  of  a  wave 
motion  whose  wave  length  becomes  Infinite.   If  7  =  c   Is  the 

o  o 

highest  critical  speed  there  Is  no  real  value  of  A,  say  A  , 
corresponding  to  a  7  value  J  =  J      such  that  7„  '^  c  .   If  there 
were  the  Inequality  (3- 17)  would  show  that  we  could  force  7  to 
zero  by  Increasing  7  to  some  value  7'.   This  would  produce  a 
critical  speed  higher  than  c   contrary  to  the  assumption  that  c 
Is  the  highest  critical  speed.   In  other  words  there  Is  no  bounded 
flow  other  than  the  equilibrium  flow  if  7  ^  c  .   It  Is  evident 
from  (3.19)  and  (3.1?)  that  the  formula 


(3.20)     gh  /2  /   ^°^^^   g  \    ^ ^dy 

_1  [c+v^(y)?  {^  P^(ti)[c  +v^{ti)]2 


provides  an  estimate  for  the  highest  critical  speed.   This 
estimate  In  general  Is  such  that  c  _^  c. 

Under  some  circumstances  the  formula  (3.20)  actually  yields 
the  highest  critical  speed.   If  the  density  Is  constant  (3.20) 
gives 


il  [c+VQ(y)]^  ^_-L  [c+v^(ti)]^ 


which,  after  an  Integration,  Is 

0 

(3.21)  gh  r  ^ ^  =  1 

•{-L  [c  +  v^(y)]^ 
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This  is  a  known  formula  for  the  critical  speeds,  c  =  c,  where 
the  density  is  constant.  For  a  discussion  of  this  formula  and 
other  ways  of  deriving  it  see,  for  example.  Burns  [4],  or 

Peters  [5].   If  we  set  v^(y)  =  0  in  (3.21)  we  find  the  well  known 

2 
result  c   =  gh  for  the  critical  speed  in  a  rectangular  channel 

when  the  density  is  constant  and  the  equilibrium  flow  is  without 

vorticity . 

If  the  density  is  not  constant  then  the  assumption  (In 

order  to  have  stability)  is  that  the  density  does  not  increase 

as  y  Increases.   Hence  we  see  from  (3.20)  that 


I.,M^/^^/  ^ ,dy 


<  gh  /2  '      ° 


-  1     /       dT) 

U   [o  +  v^(y)]^   P^^  {,  [C+V^(T1)]2  ^"^ 


1  <  gh 


0 

r  dy 

ij_    [c+v^(y)]2 


This  means  that  the  highest  critical  speed  for  the  case  of 
variable  density  cannot  be  greater  than  the  highest  critical 
speed  for  constant  density. 

When  the  flow  possesses  no  vorticity  due  to  a  velocity 
variation,  i.e.  v  (y)  =  0,  the  formula  (3.20)  reduces  to 

0      y  ^ 

t'-22)  2/  p„(y)  /  ^dy=   =    . 

_1  _1    O  g   h 
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If  the  density  variation  is  exponential,  say 


f    \  -2ky 


we  find 


-1       -1  ^  '^ 


-2k    T      h 
1_   ^  e 1  ^   c 

^   2k^    21?   g^h^ 


and  when  k  is  small 

4             2  ^ 

c  T    2k  ,  k  2k^ 

-27^-  1  -  —  +  ^  -  Y^  • 

g  h        -^     ^  -L^ 

For  k  small,  this  agrees,  up  to  and  including  second  order  terms, 
with  the  approximation  to  the  highest  critical  speed  given  in 
Peters  and  Stoker  [2],  namely 


2 

c 


1^   1  +  k/3 


In  the  latter  paper  the  critical  speeds,  c  ,  are  defined  by 

2ks 

s  -  k 
(3.23)  "^ 

2    oi   gh   ,  2 
s   =  2k  •  -^  -  k 

c 
n 


Equations  (5.23)  can  be  found  by  explicitly  solving  (3. 19) 
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It  should  be  pointed  out  that  our  results  hold  for  an 
equilibrium  flow  which  is  composed  of  homogeneous  layers.   At 
an  Interface  where  the  density  is  discontinuous,  sat  at  y  =  -r, 
the  linear  theory  requires  continuity  in  the  vertical  velocity 
component  and  the  pressure.   In  terms  of  ^(y)  the  interface 
conditions  are 

■ij/{-r  -  0)  =  i/i-r  +0) 


Po(-r-  0)[u^(-r)ty(-r  -O)  -^(-r-0)] 


=  PQ(-r+0)[u^(-r)^y(-r +0)  -  ^(-r+O)] 


It  can  be  verified  that  these  conditions  are  automatically 
satisfied  by  the  integral  equation  {5.I6).   It  is  sufficient  here 
to  confine  the  discussion  to  the  case  of  a  medium  with  just  two 
layers  and  v  =  0.   Suppose  that  the  lower  layer  is  defined  by 
-1  -  r  <  y  <  -r  where  the  density  p   is  p   =  1,  and  that  the  upper 
layer  is  defined  by  -r  <  y  <  0  where  the  density  is  p   =  p   "^  1. 
In  terms  of  the  original  variables  the  depth  of  the  lower  layer 
is  h  and  the  depth  of  the  upper  layer  is  H  =:  rh  so  that  r  is  the 
ratio  of  the  depth  of  the  upper  layer  to  that  of  the  lower  layer 
and  the  corresponding  density  ratio  is  p  /l  =  p   *=  1.   The 
integral  equation  for  V(y)  is 

0  0 

<f>{y)  =  /    k^(y,Ti;7)(t)(ii)dTi  -   ^      j         k^Cy  ,ti  ;  7)Mil  )dTl 
-1-r  -1-r 
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where 

My)  =  /pQlyJ-UQ{y)^y(y)  . 

The  equation  which  prevails  for  A  =  0  and  determines  the  critical 

speeds  Is 

0 


<t'(y)  =J        k-j^{y,Ti;7)({)(Ti)dii 


-1-r 
or 


If  V  =  0  SO  that  u  =  7//gh  we  have 

y  0 


This  leads  to  the  following.   If  -1  -  r  ^  y  "^  -r  we  have 

-r  0 

u^My)    =  J        <t>(Tl)dTi    +  Pq  J       HTl)dT] 
-1-r  -r 

and  after  Integration 

-r  -r  0 

(3.24)  ulf       <})(y)dy=jr       (l)(Tl)dTi   +  p^y^    MTl)dTi    , 

-1-r  -1-r  -r 

If   -r   <   y   <   0  we  have 
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-r 
■  1-r  -r 


-I-  u 


or 

0  -r  0 

(3.25)        v^  j     (l)(y)dy  =  rpQ    /        <t){n)dTi+r   /      <l){Ti)dTi    . 
-r  -1-r  -r 


We  cannot  have  both  of  the  integrals 

0  -r 

/  <|)(y)dy   and  J       <t)(y)dy 
-r  -1-r 

equal  to  zero  because  this  would  imply  ({)  =  0  and  t^  =  0.   Hence 
the  determinant  of  the  equations  {3. 24)  and  (3-25)  must  be  zero. 
This  gives 

(3.26)  u^-  (1+  r)u^  +  r(l-  p^)  :=  0  , 

an  equation  which  defines  two  critical  speeds.   It  is  the  same 
as  that  given  by  Peters  and  Stoker  [2].   The  higher  critical 
speed  is  given  by 


Cq        1+  r  +  /(I-  rf+  4p^ 
gh  =  2 


<   Jl  +2p^r +r^      . 


This  should  be  compared  with  the  estimate  of  the  higher  critical 
speed  which  comes  from  (3.20)  namely 
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2 

c 

gTT 


0 


y 


'f       Po(^)/   p^* 


-1-r 


-1-r 


/  -r    y  0   -r        0   y 

(  -1-r  -1-r  -r   -1-r       -r   -r 


drjdy 


l^^yrri^TTT^ 


which  confirms  the  fact  that  for  the  higher  critical  speed  c  , 

c^  Is  such  that  c   <  c. 
o  o  — 


^«   Three-dimensional  Motion.   Channel  with  Arbitrary  Cross 
Section 

This  part  of  the  paper  is  concerned  with  a  three-dimensional 
case  of  our  problem  in  which  a  liquid  of  constant  density  is  con- 
fined to  a  channel  whose  cross  section  is  like  that  shown  in 
Fig.  4.1, 

'^   s        'l' 


Figure  4.1 
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The  linearized  equations  for  the  motion  come  from  (2. 15)-(2. 20) 

by  setting 

p   =  1 
^o 


They  are 


cr  =  0 


u  = 
o 


7+  VQ(y,z) 
/gh 


(4.1) 


U  U  +  U   V  +u   w  =  -  p 
ox   oy    oz      *^x 


{ 


U  V 
O  X 


u  w 

O  X 


=  -  p. 


-  p. 


u  +  V  +  w    =0 
^    X    y    z 


At  the  channel  wall  given  by  y  =  q(z)  we  must  have 


(4.2) 


V  =  w 


dq 
dz 


If  the  equation  of  the  free  surface  Is 

y  =  f(x,z) 

the  linearized  free  surface  conditions,  to  be  satisfied  at  y  =  0, 
are 


(4.5) 


v(x,0,z)  =  u^(0,z)f^(x,z) 


-f^(x,z)  +  p^(x,0,z)  =  0 
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It  seems  that  the  easiest  way  to  conduct  an  analysis  of 
the  above  equations  Is  to  regard  the  dynamic  pressure  p  as  the 
fundamental  dependent  variable.   For  any  point  In  the  bounded 
domain  D  the  pressure  p  must  satisfy 


o  o  o 

In  addition  the  following  boundary  conditions  must  be  satisfied. 


(4.5) 


On  L:    p  =  0  . 


On  S:    Py(x,0,z)+  u^p^(x,0,z)  =  0 


We  use  V  to  denote  the  unit  outward  normal  to  the  boundary  of  D. 
If  we  attempt  to  follow  the  method  of  Welnsteln  we  are  led 
to  this  eigenvalue  problem.   For  D: 


'i)'Hi) 


o        o        o 

On  L:    f     =   0    . 

2  2 

On  S:    f     =  u  A  t//  . 
y    o   "^ 

This,  however.  Is  not  a  standard  eigenvalue  problem  because  the 

2 
eigenvalue  parameter  A  appears  In  the  boundary  condition  along  S 

(with  unusual  sign)  and  S  covers  only  part  of  the  boundary  of  D. 

Instead  of  trying  to  establish  the  existence  of  a  complete  set  of 

elgenfunctlons  for  this  case  we  will  proceed  by  using  an  alternate 
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method  based  on  the  generalized  Fourier  transform. 
Let  the  right  hand  transform  of  p(x,y,z)  be 


00 

5(A,y,z)  =Je^   p(x,y,z)dx 


where  Im  A  =  a  >  0.   Let  the  left  hand  transform  of  p  be 


u 

1  (A,y,z)  =  /  e^  ^  p(x,y,z)dx 
1         ^ 

-00 


where  Im  A  =  b  <  0.   By  taking  the  magnitudes  of  a  and  b 
sufficiently  large  these  transforms  exist  for  any  p  of  exponential 
order.   The  recovery  formula  for  p  Is 

CD  +la 
i   r       -ixA  x^^ 
■^    J  e     JdA 

-00  +la 
p(x,y,z)  = 

oo+lb 

-00  +lb 


+  ^ 


The  application  of  the  right  hand  transform  to  (4.4)  and  (4.5) 
for  X  >  0  gives 

(k   6)     ^{^]  +  S  /5y^  _  ^^  ,  Px'°;^'y'    iXptO.z.y) 


o         o      o 


to  be  satisfied  in  D,  subject  to 


(4.7)  f     =  0 

V 
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on  L,  and 

(4.8)    f  =   u^A^  +  p^(0,z,0)u^(0,z)  -  iAp{0,z,0)u^(0,z) 

on  S.   Similarly  the  application  of  the  left  hand  transform  to 

(4.4)  and  (4.5)  for  x  <  0  gives 


3i" 


/IlA  ^  S  (\^        ''%        ^^[0,^,1)        iAp(0,z,y) 
V  ^o^     ^   ""o'         % 


to  be  satisfied  in  D,  subject  to 


(b   =  0 
^IV 


on  L,  and 

¥   =  u^A^  -  p^{0,z,0)u^(0,z)  +  iAp{0,z,0)u^(0,z) 

on  S. 

From  the  equations  and  boundary  conditions  which  J  and  ^ 
must  satisfy,  it  is  evident  that 


$  {'^,z,y)   =  -  f(A,z,y) 


Therefore  we  see  that 


(4.9)         p{x,y,z)  =   -  ^     I   e"^^^  J(A,z,y)dA 


C 


where  C  is  the  path  C  =  C-]_  +C2  shown  in  Fig.  4.2.   The  lines  C^ 
and  Co  are  parallel  to  the  real  axis  in  the  A-plane  and  their 
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A-plan( 


C: 


Figure  4.2 


distances  from  the  real  axis  can  be  adjusted  to  admit  functions 
p  of  various  exponential  orders. 

An  integral  equation  formulation  for  the  determination  of 
^  can  be  used  to  show  that  it  is  expressible  as  a  ratio 

i sjpri — 

in  which  each  of  f   and  cd  is  an  entire  function  of  A.   If  the 
disposition  of  the  zeros  of  odCa)  is  known  then  the  behavior  of 
p  with  respect  to  x  can  be  found  from  (4.9)  by  using  the  theory 
of  residues.   Also,  if  we  require  p  to  be  bounded  we  must  choose 
the  path  C  in  (4.9)  so  that  it  contains  only  real  poles  of  ^/o) 
and  if  necessary  ip   must  be  modified  so  that  these  poles  are  poles 
of  the  first  order.   Now,  the  substitution  of  ^/cd  for  J  in  (4.6)- 
(4.8)  shows  that  the  zeros  of  (X)(A)  are  just  the  eigenvalues  of 
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3  /*z\    a  (*A_>^% 


o  o      o 


with 

(4.11)  f^=o 

on  L,  and 

(4.12)  ^y  =  u^A^^ 

on  S.   Hence  our  problem  Is  again  reduced  to  a  study  of  eigen- 
values.  Here,  however,  we  do  not  need  to  know  anything  about  the 
completeness  of  the  set  of  elgenf unctions. 
For  the  operator 

ais^siB 

o         o 
we  have  the  following  identities 

D  D   ^      o     -    L+S   o 

(4.14)       ff  [fEie)  -9EW]    =f     j   %  ^]ds 
D  L+S  I    ^o    ^ 

where  s  is  the  arc  length  along  the  boundary  of  D.   If  A  and  f 
are  the  respective  conjugates  of  A  and  -p   and  if  we  identify  the 
conjugates  with  d   in  the  above  identities  they  show  that 


(4.15)        (A^-A^)  j  7  fiJa.  -ji  ^  \  =  o 

D   ^o 
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and 


(4.16) 


a2 


f       #dz    -   [f    ^\=  [f 


V     ^1 


D 


D 


tl/  il/    +11/11/ 


u 


o 


From  (4.15)  and  (4.l6)  it  follows  that 


(^.17) 


(A^ 


^'l[l 


D      J- 


t1/  il/    +  il/  il/ 


u 


dzdy  =  0    , 


and  from  this  we  conclude  that  the  eigenvalues  are  either  real 
or  pure  imaginary  numbers.   It  should  be  noted  that  if  A  is  an 
eigenvalue  then  so  is  -A.   It  should  also  be  noted  that  A  =  0  is 
an  eigenvalue  of  (4. 10)-(4. 12)  and  that  the  corresponding  eigen- 
function  is  f^   =   const.  ^  0.   We  infer  from  the  last  two  observa- 
tions that  A  =  0  is  at  least  a  triple  zero  of  cd(A)  if  the  speed 
7  corresponds  to  a  transition  from  real  eigenvalues  to  pure 
imaginary  eigenvalues. 

Suppose  A  =  0  is  a  simple  zero  of  cd(A)  and  the  only  one  in 
a  strip  YZ.   which  contains  the  real  axis.   Then,  as  we  can  see 
from  (4.9),  the  only  bounded  solution  for  p  is  given  by  the 
residue  of 


e-^^^  ?(z,y;A)  = 


-ixA  , ,     , , 
e     ^(z,y;A) 


at  A  =  0;  that  is,  p  =  p(z,y).   However,  if  p  does  not  depend  on 
X  the  equations  (4.4)-(4.5)  show  that  p  is  constant.   With  p 
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constant  it  follows  from  the  equations  (4.1) -(4. 3)  that  the  flow 
must  be  an  equilibrium  flow. 

Suppose  next  that  the  speed  7  is  such  that  A  =  0  is  a  triple 
zero  of  cd(A)  and  the  only  one  in  YZ   '      ^'-'^   this  case,  the  residue 
of 

e     i{z,y;A)  =  cJ{i) 

at  A  =  0  would  generate  unbounded  terms  for  p  unless  f[z,Y;\)    =  0; 
or  unless  ^(z,y;A)  possesses  a  double  zero  at  A  =  0.   If  we  Impose 
either  of  these  conditions  we  find  again  that  p  =  p(z,y)  and  this, 
as  we  indicated  above,  implies  an  equilibrium  flow.   If  A  =  0  is 
the  only  zero  of  co(A)  in  YZ   '>    ^^^  1^  i^  ^^  ^  zero  of  odd  multi- 
plicity greater  than  three;  an  analysis  similar  to  the  above  leads 
to  the  same  result.   In  other  words,  we  conclude  from  this  para- 
graph and  the  last  one  that  the  equilibrium  flow  is  the  only 
bounded  flow  if  all  the  eigenvalues  are  pure  imaginaries  including 
A  =  0. 

With  the  real  eigenvalue  A  we  can  associate  the  real 
elgenfunction  Tp .      By  setting  0  =   iIj  In   (4.13)  we  have 


(4.18)      a:  /   t  dz  =  K   1 1     ^+11      2 


r  J  r  ^. , 

T^   u    XT      u 
q-,  D    o   D       o 


By  setting  0   =  1  in  (4.14)  we  have 
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(4.19) 


i/dz 


D    o 


An  application  of  Schwartz's  inequality  to  the  last  equation 
gives 


/  H-ff^-/f 


D    o   D 


^o 


With  this,  (4.18)  leads  to 


(4.20)  r 


D        %    q, 


t  dz  - 


r-  1' 


^  q- 


T^dz 


2  A 


^//  ^  •  ir 

D   ^o   D 


1| 


■  2    2- 

^  +  V/ 


u 


For  the  case  of  the  rectangular  channel  f   does  not  depend 
on  z  and  the  inequality  (4.20)  reads 


(4.21) 


A^^^(O) 


0 

/  % 

1  u 

-1  o 


-  1  > 


0 

-1   o 


— 2~ 

■  1    o 


This  shows,  as  we  deduced  in  a  different  way  in  Section  3>  that 
the  critical  speed  is  given  by 


'        dy 
.1  % 


and  that  no  real  non-zero  eigenvalue  exists  if  7  is  such  that 
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0 
(4.22)  f      %  1  1 

-1   o 


Hence,  as  we  have  shown  above,  for  speeds  7  which  satisfy  (4.22) 
the  only  bounded  flow  is  the  equilibrium  flow. 

When  the  cross  section  of  the  channel  is  arbitrary  we  can 
define  the  critical  speed  by  requiring 


[r 

1 

1 1 

^2 
J      t^dz    - 

r     ^2           n 
/       T^dz 

D 

"^o 

^1 

L^Ql          -■ 

-i2 


(4.23) 


This  implies  by  virtue  of  (4.20)  that 


0 


2 " 

D      o 


/[ 


or  ^  =  const.  ^0.   With  this  eigenfunctlon,  (4.23)  becomes 
(4.24)  ^[  dzdy  .  b 


D     o 


where  b  =  qp  -  q,  is  the  dimensionless  breadth  of  the  channel.   The 
critical  speed  defined  by  (4.24)  is  the  value  of  7  which  allows 
transit  from  real  eigenvalues  to  pure  imaginary  values  through 
A  =  0.   A  more  detailed  discussion  of  (4.24)  can  be  found  in 

Peters  [5]. 

o 
If  l/u  were  a  negative  parameter  then  all  of  the  eigenvalues 

would  have  to  be  pure  imaglnarles  as  (4.l6)  would  show.   Since  the 

eigenvalues  A  depend  continuously  on  l/u  we  conclude  that  when 
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l/u  Is  small  there  must  be  a  corresponding  A*  of  least  absolute 
magnitude  which  Is  either  a  pure  imaginary  or  a  real  number  which 
is  small  in  absolute  magnitude.   If  the  magnitude  of  l/u  is 
sufficiently  small  we  cannot  satisfy  (4.24)  because  then 


// 


^  <  <  <  b 


D    % 
and  therefore  A*  must  be  a  pure  Imaginary.   Now  let  7  in 


be  decreased  until  7  satisfies 

gh  ff  -i^y-^  =  b  . 

The  eigenvalue  A  must  then  pass  through  a  continuum  of  imaginary 
values  until  the  origin  is  reached.  If  we  continue  to  decrease  7 
until 

the  eigenvalue  A  passes  to  real  values. 

We  can  now  conclude  from  the  above  analysis  that  if  the 
speed  7  is  such  that 

gh  ff  _d5ay     <  b 

then  no  real  and  non-zero  eigenvalue  exists  and  consequently  the 
only  possible  bounded  flow  in  the  channel  with  arbitrary  cross 
section  is  the  equilibrium  flow  defined  in  Section  2. 
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